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ABSTRACT 



NASA has recently published data on the effects of tem- 
perature and pressure distortion on the stability of axial- 
flow compressors. Theoretical predictions of compressor 
instability based on vorticity do not explain the observed 
compressor stability of the NASA tests. Description of 
distorted flow by vorticity is extended by a simplified model 
which incorporates vortex sheets and influence coefficients. 
To determine the upstream influence of the compressor on 
flow field distortion, a. model using parallel compressor 
theory, experimental data from NASA and a pressure ratio 
equation derived by Hill and Peterson is formulated. The 
model correctly predicts compressor blade loading calculated 
by Carta. This new model accounts for upstream effects of 
the compressor and shows that the vorticity model by itself 
is incomplete. Incorporation of the new model into a dis- 
tortion index based on radial vorticity and total pressure 
change is proposed and discussed. 
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I. INTRODUCTION 



In applications of high-specific-energy turbomachines, 
such as turbo-jet engines with axial-flow compressors, it 
is necessary to operate the compressor as close to maximum 
pressure ratio as possible without inducing unstable per- 
formance known as surge. A phenomenon of axial flow com- 
pressors is sensitivity to non-uniform or distorted flow 
fields when operating near maximum pressure ratio. As air- 
craft are designed for higher speeds and more comprehensive 
operating envelopes, the possibility of flow field distor- 
tion at the compressor face is increased. At supersonic 
speeds, shock effects within the engine intake system are 
significant. At high altitudes and low Mach numbers the 
effects of low Reynolds number may cause portions of the 
flow in the compressor to revert to laminar flow. Laminar 
flow increases the tendency toward separation which accel- 
erates compressor stall. Aircraft operations at high angles 
of attack and/or yaw also can cause distortions of the flow 
field. High performance military aircraft may be subjected 
to all of the above distortions in addition to armament re- 
lated temperature distortions. Forward firing armament 
could also cause ingestion of foreign gas into the compressor 
and is discussed in Appendix A. Compressor surge is not re- 
stricted to the airborne flight envelope. Ground operations 
with some recently developed engines have resulted in 
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destructive surge while .taxiing. Figure 1 shows some of the 
general areas of the flight envelope where large stability 
margins are required. Figure 2 shows a compressor map that 
illustrates the stability margin. The associated performance - 
penalty of a large stability margin becomes more acute as 
engine performance requirements increase. 

This situation has motivated the formulation of a dis- 
tortion index that could be used in the design, development, 
and component matching of turbo-jet engines which use axial- 
flow compressors. Conceivably such an index could be used 
in engine control systems to allow operation at reduced 
stability margins without the danger of compressor stall. 

Several approaches to distortion indices have been tried. 
In fact there has been a proliferation of indices based on 
pressure distortion or temperature distortion at the face 
of the compressor. Most data collected are in the form of 
total pressure or total temperature measured at the com- 
pressor face. A more recent approach is to formulate an 
index based on vorticity. Since vorticity can be derived 
from either temperature or pressure data, a single index 
appropriate for both forms of distortion is possible. 

Farmer [Ref. 1] investigated converting total-pressure 
fluctuations at the compressor face into vorticity. Iverson 
[Ref. 2] proposed a method of transforming total-temperature 
contours into vorticity. Shoemaker [Ref. 3] expanded Farmer's 
theory and, using a computer generated map of vorticity, de- 
rived an index, based on vorticity alone, that predicted 
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Figure 1. Flight Envelopes Showing Areas of Potential 
Compressor Instability. 
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N = compressor rotation speed 

0 = actual temperature/NASA reference temperature 
6 = actual pressure/NASA reference pressure 



Figure 2. Typical Axial-Flow Compressor Performance Map. 
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stall in selected stall events. Research conducted, concur- 
rent with this paper, by Vann [Ref. 4] , refines the Shoe- 
maker vorticity index and applies it to additional stall 
events . The refined vorticity index (Vann) is an accurate 
stall predictor in 45 percent of the cases studied. 

Braithwaite, Graber and Mehalic [Ref. 5] address the 
effects of inlet temperature and pressure distortion on 
turbo-jet performance and provide experimental results with 
which the vorticity model can be compared. In the case of 
total-pressure distortion, the vorticity model is a good 
predictor of decreased compressor performance. The vorticity 
model based on total -temperature distortion does not cor- 
relate well in all cases when compared with Ref. 5. 

To gain further insight into this apparent failure of 
the vorticity model to correlate with experimental results, 
this paper will apply the vorticity model to models simu- 
lating the conditions of Ref. 5. In addition, an extended 
model that includes the engine effects on the upstream flow 
field will be introduced and a new approach to a distortion 
index will be discussed. 
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II. THE NATURE OF THE PROBLEM 



The vorticity approach of Farmer, Iverson and Fuhs 
[Ref. 6] gives the following relationships for pressure and 
temperature distortion: 



, 1 9P’ 

r U'yP'r’ ''96* 

z 



( 1 ) 



w 



t 

r 



M 9T’ 
2r ' 9l^ 



( 2 ) 



where variables with superscripts (primed) are non-dimen- 
sionalized quantities as defined below: 



r = 


.D 

2 


D = compressor diameter 


U = 


U'a 


a = average speed of sound 


CO = 


, 2a 
“ -D 


w = fluid vorticity 


Pt ■ 


p.p^ 


P^ = average stagnation pressure 






T^ = average stagnation temper- 



ature . 



The subscripts 6* and r denote circumferential and radial 
directions, respectively. A complete list of symbols is in- 
cluded on pages 8 and 9. 

From the relationships given in Ref. 6, it also can be 
shown that in the case of combined pressure and temperature 
distortions 



II II RT 

U w - U w = - ,v 

z r r z rP^ 



9Pi 

96 ’ 






19T. 



1 + 



'P r 



9T^ 



(3) 
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Order of magnitude analysis [Ref. 6] indicates that is 

of negligible size. Reference 5 compares the effects of 
pressure and temperature distortion, separately and in com- 
bination, on the performance of J85-GE-13 turbo-jet engines. 
To induce total-pressure distortions on the face of the 
compressor, screens of varying solidity and 180° circumfer- 
ential extent (semi-circle) were placed one diameter upstream 
from the compressor. The total temperature of the air was 
distorted by means of a hydrogen burner, also of 180° extent, 
located four diameters upstream from the engine inlet. The 
burners were external to the bellmouth. The test procedure 
is discussed in detail in Ref. 5. Figure 3 shows a simpli- 
fied schematic of the test set-up used. Figure 4 is a com- 
pressor map showing the degradation of performance associated 
with three combinations of pressure and temperature distor- 
tion at the engine inlet. From Fig. 4 the most adverse 
effect on compressor performance is observed to be the case 
of temperature and pressure distortions super- imposed (180° 
overlapped) . Figure 5 shoivs the circumferential profiles 
of total temperature and total pressure at the engine inlet. 

Based on the vorticity approach to distorted flow, one 
would expect that large values of would correlate with 
large loss of stall margin. Radial vorticity can be calcu- 
lated from Eqn. (3) . IVhen 9P^/3 0* and 9T^/90* have the same 
algebraic sign, one would expect the most unfavorable case. 
Reference to Fig. 5(a) shows that for the 180° over] apped 
case, 9P^/90* and 9T^/90* have opposite signs; for this case 
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Figure 3. Schematic of Test Set-Up Used to Impose Distortion Flow Field At Compressor Inlet. 
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Distortion (From Ref. 5) . 



INLET PRESSURE OR TEMPERATURE RATIO, 



o TEMPERATURE 






(C) DISTORTIONS 0 ° OVERLAPPED. 



Figure 5. Compressor Inlet Temperature and Pressure Profiles 
witli Combined 180° Extent Pressure and Temperature 
Distortions (From Ref. 5) . 
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is small, and there should be only small loss o£ stall 
margin. For the 0° overlapped case, Fig. 5(c), 3P^/36* and 
9T^/90* have the same sign, and co^ is large. The 0° over- 
lapped case should result in large loss of stall margin. 
Observing Fig. 4, it can be seen that the opposite result 
was observed experimentally. 



J 
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III. MODELS OF FLOW WITH PRESSURE AND TEMPERATURE DISTORTION 



Effects on the flow field at the face of the compressor 
when distortions are imposed upstream by screens (pressure) 
or heating (temperature) will be examined. Table 1 contains 
a list of influence coefficients for area change (dA/A) , heat 
addition (dh^/c^T) and friction (dC£/2) that will be used in 
this analysis. These influence coefficients are adapted 
from Shapiro [Ref. 7] . 



A. PRESSURE DISTORTION 

Figure 6 is a schematic of the model used for pressure 
distortion. The streamline curvature depicted is for 
graphic purposes, and it is not assumed that the streamtube 
change actually will be as drawn. The screen is assumed to 
be of 180° extent (which is consistent with the methods of 
Ref. 5) , and the flow at station 1 is assumed to be uniform. 
From continuity 



m 



A1 



= m 



A 2 



PlUfAi = mass flow rate 



or in differential form after logarithmic differentiation 




dU ^ dA 
-tJ ^ -A • 



The physical geometry requires that 

A-'A ' A-'b • 



( 4 ) 



The assumed uniform static pressure distribution requires 
that 
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Figure 6. Pressure Distortion Model (Schematic). 



AREA CHANGE 

dA 

A 



HEAT ADDITION 
d 



c T 
P 



FRICTION 
d C^ 

“2 



dU 1 

U - 1 




Y 

TW 



dp Y 

p ^ M2 ~ 1 




[ 1 + M^ (Y ■ 1) ] Y 

1 - M2 




0 




0 



TABLE 1. Influence coefficients for area change, heat addition 
and friction. 



= dP'» 

P-'A P-'B • 



(5) 



From Table 1 



dP^ _ -Y dA. 

P-'B M2 - 1 A''B 



( 6 ) 



dP. _ - Y dA. 

P-'A M^ - 1 A-'A 



[1 + m^y-i)]y 

1 - M2 2 • 

(7) 



Combining Equations (4), (5), (6), and (7) and algebraic 
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manipulation yields the following relationship for area 
change ; 



a J 



dC 



A^A 



= [1 + (Y - 1)] 



A-^B 



( 8 ) 



Using similar techniques to investigate velocity, it can be 
shown that 

+ Y ^ 

U^A M2 - 1 A^A 1 - M2 2 ’ 

and combining Equations (8) and (9) with some simplification 
reduces to 



dU. _ [mHI + Y)- ll 

U-'a L 2(1 - M2 J 2 



From Table 1 

dU. _ 1 dA. 

U-'B ■ 1 - M2 A-'B 



( 10 ) 



( 11 ) 



and the ratio of velocity changes can be expressed as 

dU, (y + 1)M^ - 1 ^^f 

U-'a _ 2 2 

dU. -1 dA. 

U-'B 1 - M2 A^B 

which, using Equation (8), reduces to 



U-'a 

MT” 

U-'B 



for Y “ 1.4. 

Equation (8) 
streamtube A will 



tion coefficient 



(Y 1)M^ - 1 _ 2.4 M^ - 1 
1 + (Y - 1)M2 1 + 0.4 M2 



f(M) 



( 12 ) 



shows that for all real valued Mach numbers, 
increase in cross sectional area (the fric- 
is always posi tive) , and streamtube B must 
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decrease a like amount as in Equation (4) . For subsonic 
flow, Ug 2 must then be greater than Equation (12) is 

finite for all real values of Mach number but changes sign 
when M = 0.6455. For values of M less than 0.6455, Equa- 
tion (12) is negative which requires the velocity change in 
streamtube, A to be negative. Conversely when M > 0.6455 , 
U ^2 ^ ^Al* Quantitatively M may be less than 0.6455 al- 
though transonic compressors have axial Mach numbers ap- 
proaching 0.65. 



B. TOTAL TEMPERATURE DISTORTION 

Figure 7 is a schematic of the total temperature distor- 
tion model. The illustrated curvature of the streamtube is 
for graphic purposes and does not represent an assumed pat- 
tern of change. At station 2 the static pressure distribu- 
tion is assumed to be uniform. This implies that station 2 
is far enough ahead of the compressor so that there is no 
upstream effect of the engine. 

Using Table 1 



dP 



Y M^ dA 



P^B 1 - M2 A^B 



) 



(13) 



^) 

P^A 



Y M^ 



M2 



A^A 



Y M' 



dh. 



M2 c T 



(14) 



Setting Equation (13) equal to Equation (14) , because of the 
uniform pressure assumption, then 

dh. 



A^B 



dA 



A^A 



c T 
P 



(15) 
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Figure 7. Total Temperature Distortion Model (Schematic) 



Since 



dA. 

A J 



A-'B 



A^A 



Equation (15) becomes 



A-^B 



A^A 



1 ‘I'-t 

ICpT 



( 16 ) 



With the addition o£ heat, decreases and increases; 

and for subsonic flow, Ug^ must be greater than Ug^. 

The effect of heat addition on the velocity is investi- 
gated similarly 

, dh^ 



dU 



U^A 



dA 



1 - M2 A'^A 



1 - M2 c T 
P 



nJ 



U^B 



1 - M 



dA^ 

2 A^B 



(17) 



(18) 



Incorporating the streamtube area changes of Equation (16) into 
Equations (17) and (18) , the velocity change is found to be 
the same for both streamtubes. 

. dh^ 



U^A 



U^B 



2(1 - M2) CpT 



(19) 



C. COMBINED TEMPERATURE AND PRESSURE DISTORTION 

Figure 8 is a schematic of the models used for combined 
pressure and temperature distortion. Cases will be considered 
for zero degrees overlap and 180° overlap. 

From Table 1 the addition of heat causes a loss of stag- 
nation pressure. 
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Figure 8. Combined Temperature and Pressure Distortion Model for 0® and ISO® Overlap. 



In the case of 0° overlap 



. _ - Y 

P. -'a " 2 



dC 



( 20 ) 



dP. 
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B 



Y M' 



dh. 



c T 
P t 



( 21 ) 



The net area change of streamtube A is, using Equations (8) 
and (16) , 



dA 



A ^ANET " M ] 



dCf ^ dh^ 

■ 2 CpT 



( 22 ) 



To obtain results for a specific case, let us assume that 
the heat and screens are such that 
dCf 1 ‘>l't 



2 c 

P t 



Consequently 



dP. 



dP. 






Pt 



) 



Pt 



Equation (22) in terms of T^ and introducing dC^/2 in ter 
of heat becomes 

dA ^ _ [1 + (y-1) M^] ^ 

A -'ANET 2 



ms 






O c T^ 
2 P 



1 



dh. 



c T, 
P t 



dA . 

A -'ANET 



1 



4 c T 
P t 



1 

2 



dC, 



( 23 ) 



For velocity change in streamtube A, 
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1 
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dU , _ 1 ' dA . . Y 

U ''a 1 - A -'ANET 1 - M2 2 



and incorporating Equation (23) simplifies to 

dU . _ 1 + 2y ‘^‘^f 

U ■ (1 - M2) 2 2 



In streamtube B, 

dU , _ -1 dA , 

U ''B 1 - M2 A ''BNET 1 - M2 c T 

P 



dU . _ 1 dA . ‘^^t 

U ''B 1 - M2 A ''ANET 1 - c T 

P 



dU , 3/2 + (y - 1)M^ ‘^^f 

U ''B 1 - M2 2 



( 24 ) 



(25) 



The ratio of velocity changes in the two streamtubes (for 
Y = 1.4) is 



U ''B _ 1.5 + 0.4 M^ 
M . 0.5 + 1.4 M^ 

U ''a 



f(M) 



(26a) 



Examination of Equation (26a) shows that the velocity in 
streamtube B is higher for all subsonic values of M. Recall 
that this result is for the specific case of equal total 
pressure in streamtubes A and B which resulted from the ar- 
bitrary relation of dC£/2 = 1/2 • dh^/CpT^ . 

For the case of 180° overlap the same method of develop- 
ment is used. Figure 8(a) is the model used, and there are 
three stations illustrated for the purpose of adding the 
solutions of each distortion. Assume that static pressure 
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is uniform at stations 2 . and 3. Adding Equations (19) and 
(10) and again considering the specific case of 



dCf 1 dht 

2 2 c 

P t 

the velocity change in streamtube A is 



U ''a 



^^f (1 + 2yM^) 

4(1 -M2) 



For streamtube B, using Table 1, Equations (8), (16), and 
(19) the velocity change is 

“U 4(1-M2) [3 + 2M^(y-1)] 



and the ratio of velocity changes is, for y = 1.4, 



U ''a _ 1 + 2 _ 1 + 2.8M2 

^ . 3+2M2 (y-1) 3 + 0.8M2 

U ''B 



f(M) 



(26b) 



Examination of Equation (26b) shows that the velocity in 
streamtube B is always higher than the velocity of stream- 
tube A. For the case of 180° overlap of distortions, the 
total pressure will not remain equal in both streamtubes 
as in the case of 0° overlap. In the case of 180° overlap, 
the change in total pressure for the distorted streamtube 
(A) will be 

-2yM=^ ^^f 
2 

for this specific case of 

dCf ^ dh^ 
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Total pressure will remain constant in the undistorted stream 
tube (B) , 

D. HEAT ADDITION AT SMALL MACH NUMBERS 

Using the temperature distortion model o£ Figure 9, the 
effects of heat addition on the flow field will be made more 
restrictive by limiting the entering Mach number to very 
small (near zero) values. This is accomplished by placing 
the burner outside the bellmouth, where velocities are small. 
This analysis also models the tests by NASA reported in 
Ref. 5. An energy balance at stations 1 and 2, assuming 
adiabatic conditions between these stations, gives 









p^Bl 2 B1 



= 



"^A*-*^p^A2 2 ^^12^ ^B*^^p^B2 2 ^B2^ 



(27) 



Note that station 1 is immediately downstream from the burner 
For streamtubes A and B then 



" I “h ■= S^A2 * I “a2 ' VtA 



and 



c T + i- U^ 
^p^Bl 2 ^B1 



r T + — U^ 
p^B2 2 ^B2 



c T^x, 
p tB 



(28) 



(29) 



Assume uniform P at stations 1 and 2, = nig and isentropic 

flow for streamtubes A and B between stations 1 and 2. With 
the previous assumption of small Mach number, conditions at 
station 1 are specified. The stagnation pressure for both 
streamtubes at either station 1 or 2 equals the ambient 
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figure 9. Temperature Distortion Model for Heat Addition at 
Low Mach Number. 
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pressure. Static pressure for both streamtubes is equal at 
station 1. Due to heat addition The following 

isentropic relationships now can be used: 






Y-1 

Y 



and 



(30) 



T = T . — 

B2 tB IP 



B2 



Y-1 

Y 



tBl 



(31) 



Using Equations (30), (31) with (28) and (29) will give the 
relationship of velocity for streamtubes A and B 






“b2 



WtA 






1 - 



1 - 



Y-1 

/p \ ^ 



Y-1 
/p \ Y 
2 1 



Since 



[p~ 



tA 



B2 tB 



P. 



P_ 

P. 



)g = f(Y,M) 



(32) 



and for hydrogen (fuel for burner in NASA test of Ref. 5) 
oxidized in excess air, - Yg> then M^2 “ ^^B2‘ 

Investigating the area change at 2 with the same assump- 
tions 
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PAU)b2 P^^)a2 



^A2 

^B2 



B2 A2 



A2 



T 



/T 



tB 



B2 /T 



tA 



»A2 



T 

^A2 ^B2 



iT^ T 
^B2 *A2 



= 1 



and using Equation (32) 



^A2 

^B2 



/T 



tA 



/T 



tB 



(34) 



Equation (34) is restrictive in tliat the assumptions limit it 
to the case of equal mass flow in streamtubes A and B. How- 
ever, parallel compressor theory uses a mass -weighted sum of 
subcompressors to determine the operating point of a distorted 
compressor, so in the case of 180° sector distortion this as- 
sumption is not unreasonable. 



36 



IV. COMPARISON OF FLOW MODELS WITH VORTICITY THEORY FOR TEM- 
PERATURE DISTORTION 



The analysis of Section III(B), total temperature dis- 
tortion, shows that the velocities of the two streamtubes 
are the same, thus making the radial vorticity zero. Set- 
ting in Equation (3) equal to zero and solving for 3P^/80* 
yields the following relationship: 



RT. 



8P. 



rP^U 39* 
t z 



c,. 9T. 

2. t 

rU 30* 
z 



1 + ^ M^^ 



Simplifying and changing T to T^ gives 



RTt dPt 



c dT^ 
p t 2 



dP. 



1 dh . c 

- izl —L -2. 

2 RT^ ^ Cp 



and since R/c^ ^Cyl^Y 



dP^ T c dh^ 

t ^ _ Yii JO t 

P. 2 R c T 



P t 



M^ = 



- yM! 



dh. 



c 'K 

P t 



(35) 



Equation (35) is equivalent to the influence coefficient of 
Table 1, thus verifying vorticity for this case. 

For the case of heat addition at low Mach number, a more 
subjective comparison must be made. Equation (3) simplified 
for temperature distortion only shows 
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( 36 ) 



c 3T^ f 

“r ■ RU^ 3^ 1 

All values in the right hand side of Equation (36) are always 
positive with the exception of the circumferential total tem- 
perature gradient. For the velocity profile shown in Figure 
9, to^ is positive when 6* = 180° and 31^/36* is also posi- 
tive. A similar comparison at 6* = 0° also shows algebraic 
agreement in that and the total temperature gradient are 
both negative . 
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V. EXTENDED MODEL TO ACCOUNT FOR COMPRESSOR INFLUENCE 



A. PARALLEL COMPRESSOR MODEL 

Reference 5 employs the parallel compressor model as a 
method of predicting the effects of distortion on stall 
pressure ratios'. Figure 10 shows a schematic of the parallel 
compressor model with some of the assumptions used. Figure 
11 shows the predicted stall points of the compressor for 
various pressure and temperature distortions using this 
model. Figures 12 and 13 show the profiles of static pres- 
sure, total pressure and total temperature at the compressor 
inlet and outlets for both pressure and temperature distor- 
tion. From Figure 12 it can be seen that static pressure 
is not uniform at the face (inlet) of the compressor. This 
was one of the assumptions used in the development of Sec- 
tion III of this paper. The development of an extended 
model in this section will consider the effects of non- 
uniform pressure profiles, 

B. PRESSURE RATIO EQUATION 

Hill and Peterson [Ref. 8] derive a compressor pressure 
ratio equation that predicts the total pressure ratio (R^) 
across an axial compressor stage as 



where the subscripts 1 and 2 refer to stations before and 



n 



c 




(37) 
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4. BOTH EXIT TO A COAmiON STATIC PRESSURE 



Figure 10. Paralle] Compressor Model (From Ref. 5) 
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LAPPED. 
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(D) P&T 90® OVERLAPPED. 



Figure 11, Application o£ Parallel Compressor Model to Predict 
Stall Points (From Ref, 5). 
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STATIC PRESSURE, TOTAL PRESSURE 
OR TEMPERATURE RATIO, 





(B) COMPRESSOR EXIT, STATION 3. 



Figure 12. Compressor Temperature and Pressure Profiles with 
180° Circumferential Temperature Distortion (From 
Ref. 5). 
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Figure 13. Compressor Temperature and Pressure Profiles Avith 
180° Circumferential Temperature Distortion fFrom 
Re f . 5 ) . 
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after the compressor stage. Figure 14 shows the relationship 
of the angles and 3^ with a compressor blade row at wheel 
speed f2r. For flow in unstalled cascades of high solidity 
(ratio of blade chord length to blade spacing in the cascade) , 
the exit angle is relatively constant with changes of 3^^. 
Blade loading is a function of, or can be correlated to, A3. 

If 32 ~ constant, then 3^^ has a significant effect on the 
stall characteristics of a cascade. Consider Equation (37) 
^^^ith a fixed 11^ and wheel speed. An increase in temperature 
must be accompanied by an increase in 3j^. An increase in 3^ 
can be obtained only by a decrease in axial velocity. This 
is shovm in Figure 15. Figure 11 shows that parallel com- 
pressor theory substantiates the decrease in axial velocity 
for this case since corrected v>/eight flow is proportional to 
axial velocity. The model of Section III.D. shows a veloc- 
ity profile reversed from the results of this section; see 
Figure 9. An assumption of Section III.D. was uniform static 
pressure at the downstream station. The model in fact pre- 
cluded any upstream effect of the compressor itself on the 
flow field. 

C. PRESSURE DISTORTION 

Figure 16 illustrates the profiles of static and total 
pressure} static and total temperature and axial velocity 
as derived from the models of Section III and as measured 
at the face of the compressor during the NASA tests of Ref. 

5. Flow properties at the constant-static-pressurc station 
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Figure 14. Cascade Geometry for Compressor Rotor. 




Figure 15. Effects of Temperature on Axial Velocity. 
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Figure 16. Derived and Measured Profiles of P, P^, U, T for Pressure Distortion. 



agree with those illustrated by Figure 6. According to 
Figure 12, the static pressure increases in the high-P^ 
streamtube and decreases in the low-P^ streamtube. The 
flow decelerates in the high-P^ streamtube and vice-versa 
for low-P^ streamtube. Also the area of the high-P^ stream- 
tube of Figure 6 must increase. 

The screen creates a vorticity sheet which is convected 
into the compressor. This sheet separates the high-P^ and 
low-P^ streamtubes. As predicted by Helmholtz's vorticity 
laws, the vortex sheet is a flow feature which persists once 
it is created [Ref. 9] . The vortex sheet decays by means 
of viscous forces. Figure 17 shows a model that uses a 
"bent vortex sheet" associated with the change of velocity 
and streamtube area for pressure distortion. Specific points 
(A, B, C, D) in relation to this vortex sheet are indicated, 
and the change in axial velocity at these points is also 
shown in Figure 17. 

Due to bending of the vortex sheet, the flow acquires a 
Uq component in the direction of rotor motion at points A 
and B and a Uq opposed to rotor motion at points C and D of 
Figure 17. This has an influence on 3^ as discussed later. 

D. TEMPERATURE DISTORTION 

Figure 18 shows the flow properties at three stations: 
immediately downstream of the burner, at a constant- static- 
pressure location within the duct and at the rotor inlet 
station. The flow properties at the constant-static-pressure 
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Figure 18. Temperature Distortion with Zero Pressure Distortion. 



station within the duct agree with Figure 9 of Section III 
and figure 13(a). 

To determine the upstream influence of the compressor, 
focus attention on Figure 11(b) which is the parallel com- 
pressor model matching Figure 18. The hot streamtube, 
identified as 1 in Figure 11(b), has sharply reduced cor- 
rected weight flow. This implies low axial velocity com- 
pared with the cold streamtube, number 2 in Figure 11(b). 
Consequently, the velocity profile reverses from fast in 
the hot streamtube at the constant-static-pressure station 
to slow in the hot streamtube at the rotor. This reversal 
is attributable to upstream influence of the compressor. 

The original vorticity sheet generated by heat addition re- 
mains. The static pressure is not uniform at the rotor, 
which was an essential element of the models of Section III. 
The effect of the compressor must be to change the areas of 
the streamtubes by an increased corrected weight flow in 
the undistorted streamtube. Figure 19 shoAvs the bent vor- 
tex sheet as applied to the temperature distortion model 
and the velocity changes at the rotor. 

E. EFFECTS OF THE BENT VORTEX SHEET ON ROTOR BLADE INCIDENCE 

ANGLE 

Carta and Adamczyk [Ref. 10] postulate that it is not 
possible to prescribe the conditions at the compressor face 
a priori on the basis of upstream distortion profiles alone. 
As discussed in Ref. 10 the upstream- conditions are strongly 
dependent on the presence of the blade row itself. Carta and 
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Figure 19. Bent Vortex Model for Heat Addition. 
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Adamczyk present a graphic comparison of angle of incidence 
distribution as it differs immediately behind the distortion 
screen (upstream) from that at the blade row inlet. This 
angle of incidence distribution of Ref. 10 is compared in 
Fig. 20 with a semi -quantitative analysis of inlet flow 
angle based on the bent vortex sheet model for pressure dis- 
tortion. This analysis is discussed in Appendix B. 

There are no measured or predicted angle of incidence 
distributions available for comparison with the temperature 
distortion model. However, the encouraging agreement in the 
case of pressure distortion prompts a similar analysis of 
temperature distortion. Figure 19 shows the change in axial 
velocity at the points A, B, B', C, D and E, and Figure 21 
shows the change in llg at these points. A semi -quantitative 
analysis similar to Appendix B gives the inlet flow angle 
(Bj^) distribution shown in Figure 22. The area above the 
nominal clean inlet value of contributes to stall. The 
large angles of B|^ in Figure 22 lead the distorted flow sec- 
tor somewhat, but the majority of the high blade loading due 
to large values of B^ lies in the distorted sector. Observ- 
ing Figure 11(b), it can be secji that the compressor with 
distorted flow stalls first, which gives qualitative verifi- 
cation of the bent vortex mode] for temperature distortion. 
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Stall Angle Parameter a Inlet Flow Angle 







V//////// ZZZZA 

180° 360' 



(a) Inlet Flow Angle Distribution vs 6’ 




(b) Stall Angle Parameter vs 0=''‘ (From Ref. 10) 

Figure 20. Effects of Pressure Distortion on Inlet Flow 
Using Bent Vortex Model. 
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Figure 21 




Tangential Velocity Distribution Due to Curva- 
ture of Vortex Sheet. 
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Inlet Flow Angle 




Figure 22. Inlet Flow Angle Distribution for Temperature 
Distortion Using Bent Vortex Model. 
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VI. DISTORTION INDEX FORMULATION 



The results o£ Section V suggest an index to predict 
compressor instability or stall that combines total pressure 
and radial vorticity. Figure 21 indicates that the area 
experiencing the highest values of inlet flow angle are those 
that also would experience positive radial vorticity. Since 
there are existing data and digital computer programs [Ref. 3] 
that compute radial vorticity as well as then it is 

easy to search the compressor face for these values. There 
are several possibilities for processing these values, once 
they are calculated, into an index. One method would be to 
compare the sum of radial vorticity and AP^ for a specified 
sector (percentage of compressor face area) with a critical 
value. Since there is no known critical value, another ap- 
proach that eventually might lead to a critical value is 
recommended. In Figure 23 the total pressures across the 
face of the compressor are shown. P^ is the average total 
pressure. If P^ - P^ is positive, then the effect is to 
increase the blade loading as shown in Figure 20. Positive 
radial vorticity also increases blade loading as shown in 
Figure 20(a), and as in Refs. 1-3 and 6. The most signifi- 
cant effect on compressor instability or stall would be the 
largest value of AP^ and radial vorticity. The pressure 
change can be normalized with P^, and an index can be formu- 
lated by summing, over a sector, the combined values of AP^ 

and 0 ) . 
r 
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Figure 23. Total Pressure Distribution at the Compressor Face. 
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where the integer k represents the number of pressure rakes 
to be searched. By rotating or stepping the sector around 
the array of probes, it is possible to compute the maximum 
valued sector at any time. By plotting the maximum value 
of the index versus time, a graph similar to Figure 24 
could be analyzed for a possible critical value. A time 
history statistical comparison of deviation from the mean 
value of the index is another possible test to be applied 
and is done by Vann [Ref. 4] after adapting this proposed 
index to existing digital computer programs and comparing 
it with other known indices. 

It is implied that the index is a positive number at 
all times. There is nothing in the formulation that re- 
quires the index to be positive, and if the sector over 
which the index is computed is small enough then it may in 
fact be negative. Throughout this paper the term stall has 
been applied to large inlet flow angles that load the blade 
row beyond the stall limit of the individual blades. If 
the inlet flow angle is too small, a phenomenon known as 
negative stall occurs. Reference 8 discusses the plienomenon 
of negative stall, but it is important to point out tliat a 
negative distortion index would not indicate a negative 
stall. In the analysis of this index by Vann [Ref. 4], 
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Figure 24. Distortion Index vs. Time. 
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a sector of 60° was used, and the index values were always 
positive. Reference 4 addresses the size of sector used to 
compute this index and other commonly used indices. 
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VII. CONCLUSIONS 



Interpretation o£ distorted inlet flow interaction with 
a compressor based on vorticity maps [Refs. 1-3] is satis- 
factory as long as the assumption of constant static pres- 
sure is valid. Experiment has shown that constant static 
pressure is not the case at the compressor face, and the 
vorticity model must be expanded to account for this. The 
failure of the vorticity model to predict the conditions 
at the face of a compressor (and thus its loss of stall mar- 
gin) in the case of total temperature distortion is a result 
of not considering the effects of the compressor on the up- 
stream flow field. The upstream influence of the compressor 
causes an area change of the undistorted streamtube that re- 
sults in greater velocity changes than those induced by the 
vortex itself. The pressure ratio equation, Equation (37), 
and parallel compressor theory adequately predict trends in 
compressor stability. 

A model based on the bending of the vortex sheet gives 
results comparable to parallel compressor theory and in the 
case of total pressure distortion qualitatively duplicates 
Carta's results. 

Carta's stall incidence angle ratio suggests a distor- 
tion index that combines radial vorticity and total pressure 
gradients since both can be related to blade loadings that 
cause stall. Such an index ]ias been developed and has been 
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tested [Ref. 4] against other distortion indices using common 
data for comparison. The new distortion index has proven 
equal to or better than the distortion indices with which it 
was compared for accuracy in predicting compressor stall. 
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APPENDIX A 



Ingestion of hot gas from forward firing armament sys- 
tems has more than temperature effects on the performance 
of an axial-flow compressor. The characteristic line of 
corrected engine rotational speed (N/ /S = constant) appears 
to be a function of temperature alone since 0 is defined 
as the ratio of actual temperature to some reference tem- 
perature, usually a sea level standard day temperature. 
However, the reference to temperature alone is a convention 
that has evolved from the fact that compressors in gas 
turbine engines for aircraft applications use only air as a 
fluid medium. The corrected engine rotational speed (N//0), 
sometimes identified as referred speed, is a similarity 
parameter designed to ensure comparability of data taken 
under varying conditions of testing. In tlie case of a com- 
pressor blade row it is desirable to preserve the blade 
angle of attack. Consider two "different” compressors as 
in Figure A-1. If the angle of attack is to remain the same 
for both compressors, then 3 must be equal to 3'. 3 is the 

inlet flow angle but differs from the angle of attack only 
by a constant related to blade geometry as discussed in Ap- 
pendix B. For similarity 



tan 3 




Q ' r' 




and M must equal M'. Therefore 
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Figure A-1. Similarity for Compressors. 



flr _ ' r ' 

aM a'M' 

where a is the speed of sound. Cancelling the Mach numbers 
and substituting N (usually RPM) for the angular speed n, 
one obtains the relationship 



= N'r' 
a a ' 



Since the compressor performance maps usually relate the 
characteristics of a specific compressor under differing 
conditions, the radius r could be eliminated and the similar 
rotational speeds expressed as 



N 

a 



N X, XT I a 

— r > or N = N ' — p 
a ' ’ a ' 



Since a^ = yRT> the ratio of sonic speeds reduces to the 
ratio of the square root of the temperature if y 3^rid R are 
the same for both compressors. This is usually the case 
since y is a constant for air within the temperature range 
encountered and y = y' is fundamental to similarity. The 
gas constant R varies with molecular weight, however; and 
if the fluid medium is air plus foreign gas (steam or pro- 
ducts of combustion for example) , then the similarity rela- 
tionship must retain sonic velocity as a variable. In fact 
some literature (in particular older literature) uses the 
corrected rotational speed notation of 

N ^ iL 
XT actual 

N = 

std a 



6.S 






% 



where L is a characteristic length, usually the tip radius 
of the first stage if comparing different sized compressors. 

All corrected rotational speeds in this paper are the 
more usual N/ /© since foreign gas distortion is not con- 
sidered. References 2 and 3 analyze the effects of foreign 
gas distortion or vorticity and how it affects compressor 
performance . 
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APPENDIX B 



Carta and Adamczyk [Ref. 10] define a stall angle 
parameter a as the ratio of instantaneous incidence angle 
to tlie stall incidence angle of the compressor blade. They 
further define the incidence angle a to be the difference 
between the stagger angle X and the inlet flow angle 3^^. 
From Figure 14 



a = 3j^ - X, or3j^ = a + X 

and 3^ differs from a only by a constant X that is a func- 
tion of blade geometry. It will not be inconsistent then 
to plot a distribution of inlet flow angle 3j^ vs instead 
of a vs 6* as in Ref. 10. 

In this semi-quantitative analysis a nominal clean in- 
let (undistorted) value of 45° will be assumed for 3^^. Tlie 
changes in inlet flow angle at the points A, B, B', C, D 
and E of Figure 17, using the change in and Ug of Figures 
17 and 21, are predicted graphically in Figure B-1. The 
predicted inlet flow angles are plotted in Figure 21 with 
the stall parameter a from Ref. 10. 

Although the magnitude of many of the variables in this 
semi-quantitative analysis are somewhat arbitrary, the di- 
rections in which they were applied have been developed 
carefully in the text of this paper. Deviations from the 
magnitude used sliould give the s.'une general shape of inlet 
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figure B-1, Compressor Effects on Inlet Flow Angle. 
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flow angle curve as that 'of Figure 21. The important point 
is that the general trend of inlet flow angle closely parallels 
that calculated by Carta and Adamczyk and reinforces the bent 
vortex model as having accounted for the upstream effect of 
the compressor on the flow field. 
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